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Abstract 

The masses of baryons containing two heavy quarks and their couphngs to the cor- 
responding quark currents are evaluated in the framework of NRQCD sum rules. The 
coulomb-like corrections in the system of doubly heavy diquark are taken into account, 
and the contribution of nonperturbative terms coming from the quark, gluon and mixed 
condensates as well as the product of quark and gluon condensates, is analyzed. The 
higher condensates destroy the factorization of baryon and diquark correlators and pro- 
vide the convergency of sum rule method. As a result the accuracy of estimates is 
improved. 
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1 Introduction 



Along with an experimental search for an explanation of electroweak symmetry breaking and 
a physics beyond the Standard Model, high energies and luminosities of particle accelerators 
running or being planned and under construction provide a possibility to observe rare pro- 
cesses with heavy quarks. A competitive topic here is a study of dynamics in flavored hadrons 
containing two heavy quarks. It can play a fundamental role in an extraction of primary 
parameters of quark interactions, since a distinction between the QCD effects inside the dou- 
bly and singly heavy hadrons allows one strictly to constrain incalculable nonperturbative 
quantities determining the isolation of pure electroweak physics. 

The real possibility of such experimental measurements was recently confirmed by CDF 
Collaboration due to the first observation of Be meson As predicted theoretically this 
long-lived state of b and c quarks has the production cross sections, mass and decay rates, which 
represent characteristic values for the doubly heavy hadrons. Thus, the experimental search 
for the doubly heavy baryons can also be successful. Of course, such the search would be more 
strongly motivated if it would be supported by modern theoretical studies and evaluations of 
basic characteristics for the doubly heavy baryons. 

Some steps forward this program were already done. First, the production cross sections 
of doubly heavy baryons in hadron collisions at high energies of colliders and in fixed target 
experiments were calculated in the framework of perturbative QCD for the hard processes and 
factorization of soft term related to the nonperturbative binding the heavy quarks . Second, 
the lifetimes and branching fractions of basic decay modes were evaluated in the Operator 
Product Expansion combined with the effective theory of heavy quarks, which results in series 
over the inverse heavy quark masses and relative velocities of heavy quarks inside the doubly 
heavy diquark . Third, the families of doubly heavy baryons, which contain a set of narrow 
excited levels in addition to the basic state, were described in the framework of potential 
models so that the picture of spectra is very similar to that of heavy quarkonia. Fourth, 
the QCD sum rules |^ were explored for the two-point baryonic currents in order to calculate 
the masses and couplings of doubly heavy baryons 0. However, the latter analysis contains 
some disadvantages, which are related to instability of sum rules in a region of parameters 
defining the baryonic currents. This fact leads to quite large uncertainties in the calculations. 

In the present paper we analyze the NRQCD sum rules for the two-point correlators of 
currents corresponding to the doubly heavy baryons. The basic physical motivation of such 
consideration is a nonrelativistic motion of two heavy quarks inside a small size diquark bound 
with a light quark. This fact leads to the definite expressions for the structure of baryonic 
currents written down in terms of nonrelativistic heavy quarks. To the leading order of in- 
verse heavy quark mass and relative velocity of heavy quarks inside the diquark, the NRQCD 
sum rules require the account for hard gluon corrections to relate the nonrelativistic heavy 
quark correlators to the full QCD ones. The corresponding anomalous dimensions of baryonic 
currents were calculated with the two- loop accuracy 0. The NRQCD structure of currents 
corresponds to a definite choice of parameters in the full QCD expressions. Those values of 
parameters are inside the instability region, observed in the analysis done previously 0. We 
find the simple physical reason for this loose of convergency: the behaviour of quantities versus 
the sum rule parameters (the Borel variable or the moment number of spectral density) is de- 
termined by the presence of doubly heavy diquark inside the baryon and, hence, the difference 
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between the masses of baryon and diquark. This mass difference takes a dominant role unless 
we introduce the corrections corresponding to the essential nonperturbative interactions be- 
tween the doubly heavy diquark and the light quark composing the baryon. In NRQCD sum 
rules this introduction is realized in terms of nonperturbative condensates caused by higher 
dimension operators. We show that the stability of sum rules can be reached due to the 
account for the product of quark and gluon condensates in addition to the quark, gluon and 
mixed condensates. This product was not taken into account in the previous analysis in full 
QCD. Moreover, we carefully take into account the coulomb-like Usjv corrections inside the 
heavy diquark, which enhances the relative weight of perturbative parts with respect to the 
condensate terms in the calculated correlators. 

In Section 2 we define the currents and represent the spectral densities in the NRQCD 
sum rules for various operators included into the consideration. Section 3 is devoted to the 
numerical estimates. We find the masses of basic states, which are close to the values obtained 
in the potential models. The results are summarized in Conclusion. 



2 NRQCD sum rules for doubly heavy baryons 
2.1 Baryonic Currents 

The currents of baryons with two heavy quarks S*^, and H'^J^, where o denotes various 
electric charges depending on the flavour of light quark, are associated with the spin-parity 
quantum numbers = l"*" and = for the heavy diquark system with the symmetric 
and antisymmetric flavor structure, respectively. Adding the light quark to the heavy quark 
system, one obtains = for the E'^^ baryons and the pair of degenerate states = 
and = I"*" for the baryons S^^, S^^, S^^ and S*^, S^*, S^^. The structure of baryonic currents 
with two heavy quarks is generally chosen as 

J = [Q'^CTrQYq'e,,^. (1) 

Here T means transposition, C is the charge conjugation matrix with the properties CjJ^C'^ = 
—7^ and CjJC'^ = 75, i,j,k are colour indices and r is a matrix in the flavor space. The 
effective static field of the heavy quark is denoted hj Q. To the leading order over both 
the relative velocity of heavy quarks and their inverse masses, this field contains the "large" 
component only in the hadron rest frame. 

Here, unlike the case of baryons with a single heavy quark ^ , there is the only independent 
current component J for each of the ground state baryon currents. They equal 

J^t,^ = [Q'''CT^"'Q']-Jml6q%k, (2) 

JS^o^ = [Q'''CTrQ']q%k + lr[Q''^C^"'Q']-w%k, 

where JS*^^ satisfies the spin-3/2 condition 7„J5«>^ = 0. The flavor matrix r is antisymmetric 
for SJ,* and symmetric for Sgg and Hqq. The currents written down in Eq. are taken in 
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the rest frame of hadrons. The corresponding expressions in a general frame moving with a 
velocity four- vector can be obtained by the substitution of 7"^ — ;> = — -^v^. 

To compare with the full QCD analysis we represent the expression for the Jg/o current 
given in 0] 

so that the NRQCD structure can be obtained by the choice of ri = r2 = 1 and = Q 
and the antisymmetric permutation of c and b flavors. As we have already mentioned in 
the Introduction, the authors of [|^ reported that the convergency of OPE is "bad" in the 
region around the point given by the NRQCD parameters. The instability results in huge 
uncertainties of estimates. To get rid this disadvantage, we analyze the NRQCD sum rules in 
details. 



2.2 Description of the method 

In this section we describe steps required for the evaluation of two-point correlation functions 
in the NRQCD approximation and the connection to physical characteristics of doubly heavy 
baryons. We start from the correlator of two baryonic currents with the half spin 

n(t) = I j dWp^(o|r{j(x), J(o)}|o) = i)Fi{t) + F2(t), (3) 

where t = k ■ v, kfj_ and denote the residual and full momenta of doubly heavy baryon, 
respectively, is its four-velocity, which are related by the following formula: 

PtM = kf, + (mi + ms)^^, (4) 

where mi and m2 are the heavy quark masses. The appropriate definitions of scalar formfactors 
for the 3/2-spin baryon are given by the following: 

n^uit) = I j rfWp^(o|r{ J.(o)}|o) = -g^i^hit) + hit)] + . . . , (5) 

where we will not concern for contributions with distinct lorentz structures. The scalar cor- 
relators F can be evaluated in deep euchdean region by employing the Operator Product 
Expansion (OPE) in the framework of NRQCD for the time ordered product of currents in 
Eqs.(|), (i, say, 

Fr,2{t) = Y,d^''\t)0,, (6) 

d 

where Od denotes the local operator with a given dimension d: Oq = 1, O3 = (gg), O4 = 
(^G^), . . . , and the functions Cd{t) are the corresponding Wilson coefficients of OPE. In this 
work we include the nonperturbative contributions given by quark, gluon and mixed conden- 
sates. So, evaluating the contribution of quark condensate operator, we use the following OPE 
for the correlator of two quark fields |T0|: 



m{qt{x)q){m) = -^^^%{^^)-[^ + ^ + ^{^G') + ...], (7) 
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where the value of mixed condensate is parameterized by introducing the variable tUq, which 
is numerically determined as ttiq ^0.8 GeV^. 

For the sake of simplicity, we write down the Wilson coefficients of unity and quark-gluon 
operators by making use of dispersion relations over t, 



TC Jq UJ — t 

where p denotes the imaginary part in the physical region of NRQCD. Thus, the calculation 
of Wilson coefficients of operators under consideration is transformed into the evaluating the 
corresponding spectral densities. 

To relate the NRQCD correlators to the real hadrons, we use the dispersion representation 
for the two point function with the physical spectral density, given by the appropriate reso- 
nance and continuum part. The coupling constants of baryons are defined by the following 
expressions: 

(0|J(x)|H^qq(p)) = ^ZHo^«(t;,MH)e^^^ 
(0|J"(x)|S-q(p,A)) = zZH.o^n™(t;,MH)e^^^ 

where the spinor field with four- velocity v and mass Ms satisfies the equation '}/)u{v^M-e) = 
u{v, M=), and u'^{v, M=) denotes the transversal spinor, so that (7™ — v'^i/})u™'{v, Ms) = 0. 

We suppose that the continuum density, starting from the threshold oJcont is equal to that 
of calculated in the framework of NRQCD. Then in sum rules equalizing the correlators, 
calculated in NRQCD and given by the physical states, the integrations above oJcont cancel 
each other in two sides of relation. This fact leads to the dependence of calculated masses and 
couplings on the value of oJcond- Further, we write down the correlators at the deep under- 
threshold point of to = + ^2) + ^ with t ^ 0, which corresponds to the limit of ^ 0. 

The approximation of single bound state leads to the following expression for the resonance 
contribution: 

which can be expanded in series Q over t. Then, the sum rules state the following equalities 
for the terms standing in front of powers of t 

[uj + mi + m2) 



TT 







where pj contains the contributions given by various operators in OPE for the correspond- 
ing scalar correlators Fj. Introducing the following notation for n-th moment of two point 
correlation function 

1 r'^""' p{uj)du 

^To the nonrelativistic approximation wc have to substitute for ^^/_^2 by the single pole in the physical 
region over t, pgr^tjpg^rff) ~ M^{M^-t) ■ 



M. = -J v... . (11) 
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we can obtain the estimates of baryon mass M=, for example, as the following: 

MbN = (12) 

and the coupling is determined by the expression 

\Zs[n]\' = MnMS^\ (13) 

where we see the dependence of sum rule results on the scheme parameter. Therefore, we tend 
to find the region of parameter values, where, first, the result is stable under the variation of 
n, and, second, the both correlators Fi and F2 reproduce equal values of physical quantities: 
masses and couplings. The problem of full QCD was the second point: the difference between 
the evaluated variables. 

The similar procedure can be described in the Borel scheme, wherein the consideration of 
Ain is replaced by the analysis of function B[w], which is defined by 

B[w] = - / p{uj)duje-^^+"''-^"'^^/'^, (14) 



equal to the resonance term given by 

Be[w] = iZspe-^^^/'^. (15) 

Then, again the scheme dependence of masses and couplings appears because of variation 
versus the Borel parameter w. 

2.3 Calculating the spectral densities 

In this subsection we present analytical expressions for the perturbative spectral functions 
in the NRQCD approximation. The evaluation of spectral densities involves the standard 
use of Cutkosky rules |TT| with some modifications motivated by NRQCD. We explore the 



prescription that the discontinuity of two-point functions under consideration can be evaluated 
using the following substitutions for heavy and light quark propagators, correspondingly: 

1 

heavy quark: ^ 2m ■ 6{po - {m + —)), 

light quark : — 27X1 ■ 5{p^ — m^). 

For the perturbative spectral densities pi^h{uj) and P2,h{^^) standing in front of unity operator 
in Fi and F2, respectively, where H = Sqq,, Sqq or Hgg, we have the following expressions: 

pi=.. (iu) = — ^ -1672^ ""^""^ y^\''/^[A29ml + A29ml + 

'^W 150157r3(u; + mi + m2)3 \mi + 7112 J 

715m^tu + 403m2tu2 + 77uj^ + 143m^(9m2 + 5tu) + 13mi(99m^ + (16) 

110m2tu + 33cj2)], 

Pi=-^^(^) = 3pi,s-^,(^) = 3pi,H- (a;), (17) 

P2-J^) = P2,H- (^)=P2,H- (^)=0. (18) 
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Note that in the leading order of theory with the effective heavy quarks their spins are de- 
coupled from the interaction, that causes the spin symmetry relations given in Eg. ([T7|) . The 
factor 3 stands because of the normalization of vector diquark current. We see also that in 
the leading approximation of perturbative NRQCD the F2 correlators are equal to zero. This 
is because the interaction of massless light quark with the doubly heavy diquark is switched 
off in this order, and there is no mass term structure in the correlator. 

The coulomb-like interaction inside the heavy diquark can be taken into account by the 
introduction of Sommerfeld factor C for the diquark spectral densities before the integration 
over the diquark invariant mass to obtain the baryon spectral densities, so that 

Pdiquark = Pdiquark ' ^5 (19) 

with 

^ 2TTas r 

— exp 



3v 
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2'Kas 

3Vi2 



(20) 



where denotes the relative velocity of heavy quarks inside the diquark, and we have taken 
into account the color anti-triplet structure of diquark. The relative velocity is given by the 
following expression: 



^12 = a/1- ^2 J (21) 



4:17111712 

- (mi - m2) 

where is the square of heavy diquark four-momentum. In NRQCD we take the limit of low 
velocities, so that denoting the diquark invariant mass squared by = (mi -|- m2 -f- e)^ and 
the reduced quark pair mass by mi2 = mim2/(mi + m2), we find 

C - ^"""^ _ e 

t>Vvi 2mi2 

at e <^ mi2. The corrected spectral densities are equal to 

pf^„ iu) = — ( ^^^^ ^ cu^[2ml + 2m^ + cul^ 

(uj) = SpC^^ (^) = /^x /22) 

Further, the spectral functions, connected to the condensates of light quarks and gluons, 
can be derived. For the quark condensate term we have the following expressions: 

P2H'o ,M = V^, 23 

' QQ' IT \1Tli + m2 J 

pl%. (u) = 3pl%,. ^{u;) = 3pl% . {u), (24) 

' QQ QQ' QQ 

' QQ QQ' ' QQ 

which have to be multiplied by the Sommerfeld factor C, wherein the variable e is substituted 
by u, since in this case we have no integration over the quark-diquark invariant mass. 

It is interesting to stress that in NRQCD the light quark condensate contributes to the 
F2 correlators, only. This fact has a simple physical explanation: to the leading order the 
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light quark operator can be factorized in the expression for the correlator of baryonic currents. 
Indeed, we can write down for the condensate contribution 

(0|r{ J(x), J(0)}|0) = (0|TK(x)gf (0)|0) • 1 ■ (0|r{ Ji(x), Ji(0)}|0) + . . . , 

where J;^(a;) denotes the appropriate di quark current with the color index j, as it is defined by 
the baryon structure in Eqs. (0). So, we see that the restriction by the first term independent 
of X in the expansion for the quark correlator in results in the independent contribution of 
diquark correlator to the baryonic one. Then, since the diquark correlator is isolated in F2 from 
the baryonic formfactor Fi, the NRQCD sum rules lead to the evaluation of diquark mass and 
couplings from F2, and estimation of baryon masses and couplings from Fi. These masses and 
couplings are different. The positive point is the possibility to calculate the binding energy for 
the doubly heavy baryons A = Ms — Mdiquark- The disadvantage is the instability of NRQCD 
sum rules at this stage, since the various formfactors or correlators lead to the different results. 
In sum rules of full QCD various choices of parameters in the definitions of baryonic currents 
result in an admixture of pure diquark correlator in various formfactors, so that the estimations 
acquire huge uncertainties. Say, the characteristic ambiguity in the evaluation of baryon mass 
in full QCD is about 300 MeV, i.e. the value close to the expected estimate of A. The analysis 
in the framework of NRQCD makes this result to be not unexpectable. Moreover, it is quite 
evident that the introduction of interactions between the light quark and the doubly heavy 
diquark destroys the factorization of diquark correlator. Indeed, we see that due to the higher 
terms in expansion (0), the diquark factorization is explicitly broken, which has to result in the 
convergency of estimates obtained from Fi and F2. Below we show numerically that this fact 
is valid. Technically, we point out that the contribution to the moments of spectral density, 
determined by the light quark condensate including the mixed condensate and the product of 
quark and gluon condensates, can be calculated after the exploration of (|^), so that 



M 
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{n + 2)! ml ^ {n + A)\ tt^ ^g,^,^' 



n\ 16 n\ 288 ' tt 

For the corrections due to the gluon condensate we have 



M^^'^l (26) 



m\m2 , , 



p^„,<> (uj) = 5^ (28mf + 41mim2 + 28mn) ■ 

'^^'"QQ'' 1344y27r(mi + m2)2(cu + mi +m2)3 

(tu^ — 7ml + 7mluj — 21m\m2 + lm2uP' + Tmao; — lm?2 + Imxu? (27) 

-^\Am\m2<jJ — 21mim2), 

pfL (uj) = 3pfL (uj) = 3pfl*o (uj), (28) 

p^l, {uo) = pf^,. (cu) = pf^.o (cu) = 0, (29) 

which are written down for O4, having the form O4 = (— G^). 

For the product of condensates wherein the gluon fields are connected to the 

heavy quarks in contrast to the light quark, we have computed the contribution to the two- 
point correlation function itself. It has the following form: 



Ffi;, {uj) = ^- , 30 



8 



H = 3F|g,\M = 3F|f : M, (31) 

' QQ QQ' 

Fl€ {uj) = Ff'S (cu) = Fl€. {uj) = 0, (32) 

where lo = —{mi + 1712) + 1 at the point under consideration, and the correlators have to be 
expanded in series of t, which gives the moments. 

Thus, we provide the NRQCD sum rules, where we take into account the perturbative 
terms and the vacuum expectations of quark-gluon operators up to the contributions by the 
light quark condensate, gluon condensate, their product and the mixed condensate. Note, 
that the product of condensates is essential for the doubly heavy baryons, and we present 
the full NRQCD expression for this term, including the interaction of nonperturbative gluons 
with both the light and heavy quarks. The correct introduction of coulomb-like interactions 
is done for the perturbative spectral densities of heavy diquark, which is important for the 
nonrelativistic heavy quarks. Finally, we find the spin-symmetry relation for the baryon 
couplings in NRQCD 

|ZhP = 3|Zh,P = 3|Zh.|^ 



2.4 Anomalous dimensions for the baryonic currents 

To connect the NRQCD sum rules to the quantities in full QCD we have to take into account 
the anomalous dimensions of effective baryonic currents with the nonrelativistic quarks. They 
determine the factors, which have to multiply the NRQCD correlators to obtain the values in 
full QCD. Indeed, to the leading order of NRQCD we have the relation 

where the coefficient Cj depends on the normalization scale /i and obeys the matching condi- 
tion at the starting point of /xq = fni + 7712- The anomalous dimensions of NRQCD currents 
are independent of the diquark spin structure in the leading order. They are equal to M 



_ dlnCj{as,f7) fas 



m 



(iln(u) \47r/ 

\/ / m=l 

7^^) = ( -2CB{3a-3) + 3CF{a-2)), (33) 



7(2) = i(-48(-2 + 6C(2))C| + C^((104- 240C(2))CB-101C,.) 
o 

-GACbNfTf + Cf{-9Cf + 52NfTf)), 

where Cf = {N^ - l)/2N,, Ca = N,, Cb = {N, + 1)/2N,, and Tf = 1/2 for = 3, A^^ being 
the number of light quarks. In Eq. (|33D we give the one-loop result with the arbitrary gauge 
parameter a, and the two-loop anomalous dimension is represented in the Feynman gauge 
a = 1. So, numerically at Nf = 3 and a = 1 we find 

7^^) = -4, 7(2) ^ -188.24. (34) 

In the leading logarithmic approximation and to the one-loop accuracy, the coefficient Cj is 
given by the expression 

CA.^-{^f (35) 
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where Po — llNc/3 — 2Nf/S = 9. To evaluate the two- loop expression for Cj we have to know 
subleading corrections in the first ag order in addition to the anomalous dimensions. These 
corrections are not available yet, so we restrict ourselves by the one- loop accuracy. 

Further, we have to determine the normalization point for the NRQCD estimates ^ = 
Hi. We put it to the average momentum transfer inside the doubly heavy diquark, so that 
/if — Trf2mi2, where denotes the kinetic energy in the system of two heavy quarks, which 
is phenomenologically independent of the quark flavors and approximately equal to 0.2 GeV. 
Then, the coefficients Cj are equal to 

Cs,,«^1.95, Ch,, ^1.52, ^ 1.30, (36) 

with the characteristic uncertainty about 10% because of the variation of initial and final 
points /io,i- 

Finally, we emphasize that the values of Cj do not change the estimates of baryon masses 
calculated in the sum rules of NRQCD. However, they are essential in the evaluation of baryon 
couplings, which acquire these multiplicative factors. 

3 Numerical estimates 

In the present paper we explore the sum rule scheme of moments for the spectral functions of 
correlators. We have to stress that in this scheme the dominant uncertainty in the results is 
caused by the variation of heavy quark masses. In the analysis we chose the following region 
of mass values: 

mfc = 4.6 - 4.7 GeV, = 1.35 - 1.40 GeV, (37) 

which is ordinary used in the sum rule estimates for the heavy quarkonia. Next critical 
point is the value of QCD coupling constant determining the coulomb-like interactions inside 
the doubly heavy diquark. Indeed, it stands linearly in front of the perturbative functions 
of diquark contributions. Thus, the introduction of cts/v-corrections is essential for both 
the baryon couplings and the relative contributions of perturbative terms and condensates 
to the baryon masses. To decrease the uncertainty we impose the same approach to the 
heavy quarkonia, where it is well justified, and then, we extract the characteristic values 
for the heavy-heavy systems from the comparison of calculations with the current data on 
the leptonic constants of heavy quarkonia, which are known experimentally for cc and bb or 
evaluated in various approaches for be. So, our calculations give the following couplings of 
coulomb interactions 

as{bb) = 0.37, q;,(c6) = 0.45, a,(cc) = 0.60. (38) 

The dependence of estimates on the value of thershold for the continuum contribution is not 
so valuable as on the quark masses. We fix the region of ujcont as 

u^t = 1-3 - 1.4 GeV. (39) 
For the condensates of quarks and gluons the following regions are under consideration: 
(gg) = -(250-270 MeV)^ = 0.75 - 0.85 GeV^ (— G^) = (1.5 - 2) • 10"^ GeV^ (40) 
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So, we have described the choices of parameters. 

Figs. represent the calculated difference of masses extracted from the Fi and F2 
correlators^ for the baryons S^c, "^bc and S^fe. We certainly see that at low numbers of moments 
for the spectral densities, the baryon-diquark mass difference can be evaluated as 

A = 0.40 ± 0.03 GeV, (41) 

which is quite a reasonable value, being in a good agreement with the estimates in the heavy- 
light mesons. In the region of mass difference stability we can fix the number of moment for 
the spectral density, say, n = 27±1 for Hftc, and calculate the corresponding masses of baryons, 
which are equal to 

Mb,, = 3.47 ± 0.05 GeV, Mh,^ = 6.80 ± 0.05 GeV, Ms,, = 10.07 ± 0.09 GeV, (42) 

where we do not take into account the spin-dependent splitting caused by the Os-corrections to 
the heavy-light interactions, which are not available yet. The uncertainties in the mass values 
are basically given by the variation of heavy quark masses. The convergency of NRQCD sum 
rules allows one to improve the accuracy of estimates in comparison with the previous analysis 
in full QCD 0. The obtained values are in agreement with the calculations in the framework 
of nonrelativistic potential models 0. 




Figure 1: The difference between the Scc-baryon masses calculated in the NRQCD sum rules 
for the formfactors Fi and F2 in the scheme of moments for the spectral densities. 



^In these figures we have fixed the value of gluon condensate (^G^) ~ 1.7 • 10^^ GeV* and arranged TOq 
in the above region to reach zero differences between the masses, though the variation of parameters leads to 
errors in the estimates quoted below. 
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Figure 2: The difference between the S^c-baryon masses calculated in the NRQCD sum rules 
for the formfactors Fi and F2 in the scheme of moments for the spectral densities. 




Figure 3: The difference between the Sft^-baryon masses calculated in the NRQCD sum rules 
for the formfactors Fi and F2 in the scheme of moments for the spectral densities. 
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Figs. show the dependence of baryon couphngs calculated in the moment scheme of 
NRQCD sum rules. We find that the stability regions for |^hP determined from the Fi and F2 
correlators coincide with those of the mass differences. So, the baryon couplings in NRQCD 
are equal to 

(1.7 ±0.3) • 10"^ GeV^ 

(3.7±0.5) ■ 10"^ GeV^ (43) 
(1.5 ±0.3) ■ 10-2 GeV^ 

The values given above have to be multiplied by the Wilson coefficients coming from the 
expansion of full QCD operators in terms of NRQCD fields, as they have been estimated by 
use of corresponding anomalous dimensions. This procedure results in 

(6.5 ± 1.2) • 10'^ GeV^ 

(8.5 ± 0.9) • 10"^ GeV^ (44) 
(2.5 ± 0.3) ■ 10-^ GeV^ 

which are inside the regions given in the analysis of sum rules in full QCD. 
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Figure 4: The couplings \Z^^'^^\^ of Scc-baryon calculated in the NRQCD sum rules for the 
formfactors Fi and F2 (solid and dashed lines, correspondingly) in the scheme of moments for 
the spectral densities. 

For the sake of comparison, we derive the relation between the baryon coupling and the 
wave function of doubly heavy baryon evaluated in the framework of potential model in the 
approximation of quark-diquark factorization. So, we find 

|ZP^|=2v^|vl>,(0)-vl>,(0)|, (45) 



7NR12 



NR|2 
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Figure 5: The couplings of H^c-baryon calculated in the NRQCD sum rules for the 

formfactors Fi and F2 (solid and dashed lines, correspondingly) in the scheme of moments for 
the spectral densities. 
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Figure 6: The couplings of Sftfe-baryon calculated in the NRQCD sum rules for the 

formfactors Fi and F2 (solid and dashed lines, correspondingly) in the scheme of moments for 
the spectral densities. 
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where \&d(0) and \l//(0) denote the wave functions at the origin for the doubly heavy diquark 
and hght quark-diquark systems, respectively. In the approximation used, the values of ^'(0) 
were calculated in B in the potential by Buchmiiller-Tye [O], so that 



4^1^^(0)1 = 0.53GeV^/^ 

4^|^^^(0)| = 0.53GeV^/^ 

4^|^j,^(0)| = 0.73 GeV=^/^ 

4^|^^j,(0)| = 1.35 GeV^/^ 



These parameters result in the estimates in the static limit of potential models 

iZf^P = 6.0 ■ 10-3 GeV*^ 



\Z^y = l.l•10-'GeV^ (46) 
= 3.9- 10-2 GeV^ 



'bb 



As it is well known from the analysis of heavy-light mesons, the approximation of potential 
models for the static limit of coupling results in the overestimation of leptonic constants of 
heavy mesons with the single heavy quark because of the large corrections in the expansion of 
corresponding currents through the effective static fields. So, the contribution by the higher 
dimensional operators is essential, and it leads to the suppression factor of about 1/2 for the 
couplings. The introduction of this factor allows us to calculate the leptonic constants of 
heavy mesons with the same wave function of light quark 

fo = 185 MeV, Jb = 115 MeV, 

which are quite reasonable estimates, being in agreement with the results of QCD sum rules 



T3[| . Then, we use the same factor for the renormalization of baryon couplings and find 

(47) 



^■PM 2 


= 1.5 


10" 


-3 


GeV 


1 ^bc 1 


= 2.8 


10" 


-3 


GeV 


1 ^bb 1 


= 9.8 


10" 


-3 


GeV 



Comparing the values in (^) with ones in (^), we see a good agreement of NRQCD results on 
the baryon couplings with estimates of potential models improved by the correction factor, if 
we take into account the ordinary accuracy about 30% for the phenomenological quark models. 

Finally, we suppose that the corrections coming from the higher orders of NRQCD expan- 
sion to the baryon couplings are not greater than 15%, since the diquark masses are quite 
large. 

Thus, we obtain the reliable estimates of masses and couplings for the doubly heavy baryons 
in the framework of NRQCD sum rules. 

4 Conclusion 

We have considered the NRQCD sum rules for the two-point correlators of baryonic currents 
with two heavy quarks. The nonrelativistic approximation for the heavy quark fields allows one 
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to fix the structure of currents and to take into account tlie coulomb-like interactions inside 
the doubly heavy diquark. Moreover, we introduce the higher order operators responsible 
for the quark-gluon condensates to get the convergency of sum rule method for two scalar 
functions of correlators. To the leading approximation, including the perturbative term and 
the contributions of quark and gluon condensates, the correlators of three-quarks state and 
the doubly heavy diquark are factorized in separate functions, so that the sum rules result in 
the different values of masses and couplings. This fact indicates the divergency of approach 
unless the product of quark and gluon condensates and the mixed condensate are taken into 
account. Then, the interaction of two heavy quarks and light quark destroy the factorization, 
which allows one to get meaningful estimates of masses and couplings. Moreover, we have also 
calculated the binding energy of doubly heavy diquark, which is in a good agreement with the 
estimates in the framework of potential models. 

Thus, the NRQCD sum rules allow us to improve the analysis of masses and couplings for 
the doubly heavy baryons and to obtain reliable results. 
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Center of Fundamental Physics in Moscow, International Science Foundation, and INTAS- 
RFBR-95I1300 grants. 
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